Dynamic Control of Collapse in a Vortex Airy Beam 
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We study the self- focusing dynamics and collapse of vortex Airy optical beams in a Kerr medium. 
The collapse is suppressed compared to a non- vortex Airy beam in a Kerr medium as a result 
of the existence of vortex fields. The locations of collapse depend sensitively on the initial power, 
vortex order, and modulation parameters. Unlike the collapses reported before for any beam, the 
collapse may occur in a position where the initial field is nearly zero while no collapse appears in 
the region where the initial field is mainly distributed. This study sheds light on how to control 
and manipulate the location of collapse based on the initial power, vortex order and modulation 
parameter. 

PACS numbers: 42.65.Jx; 42.65.Sf 

Nonlinear wave collapse has been investigated in many areas of physics, including optics, fiuidics, plasma physics, 
and Bose-Einstein condensates [1]. In nonlinear optics, collapses of beams with various spatial distributions have been 
studied [1-3]. One of the challenges facing this field of interest is the possibility of controlling and manipulating the 
collapse dynamics. Recently, the Airy beam has attracted considerable attention after its experimental generation 
by Siviloglou et al. [4] due to its intriguing properties and potential applications such as weak-diffraction, transverse 
acceleration [4-6], self-healing [7], and sorting microscopic particles [8]. The evolution characteristics of an Airy 
beam in nonlinear medium have been studied [9-17], such as plasma channel generation [13], laser filamentation [14], 
supercontinuum and solitary wave generation [15-17]. Vortices have been the subject of many studies and appear in 
many branches of physics [18]. A vortex Airy beam is formed by superposition of an Airy beam and a vortex optical 
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field. Recently, interesting propagation dynamics and non-classicality of a vortex Airy beam have been reported [19, 
20]. The nonlinear dynamics of a vortex Airy beam in a Kerr medium is also expected to give some interesting and 
nontrivial properties. 

In this letter, we investigate the spatial collapse dynamics of a vortex Airy beam. The coupling between the vortex 
and Airy beam strongly affects the nonlinear dynamic properties of the vortex and Airy beam in the Kerr medium. 
The diffraction against self-focusing of a vortex Airy beam is effectively enhanced compared to a non-vortex Airy 
beam. In addition, the vortex of beam tends to suppress the collapse in a Kerr medium. For a non- vortex Airy beam 
in a Kerr medium, partial collapse can occur in the beam's center [12]. However, the collapse of the vortex Airy beam 
never occurs at the beam's center. This is because the center of the vortex is located at the center of the beam. Our 
findings reveal that the position of the partial collapse and the propagation distance for the appearance of partial 
collapse are dependent on the initial powers, vortex orders, and modulation parameters. The partial collapse may 
occur in the main lobe, side lobes, or outermost lobes, depending on the strength of initial powers, the order of vortex 
and the modulation parameters. The collapse can occur in the position where the initial field is almost zero while 
no collapse appears in the position where the field originally exists with a appropriate power. We further show that 
the initial power, vortex order and modulation parameters can be exploited to control and manipulate the position of 
collapse dynamics in a vortex Airy beam. Even if the initial power is ultrahigh, the partial collapse occurs separately 
at the side lobes and the beam still propagates along a similar accelerating curve trajectory as that of an Airy beam 
in free space. Since the field distribution of a vortex Airy beam is modulated by some exponential factors, we also 
study the effect of these exponential factors on the nonlinear evolution of the beam in a Kerr medium. Finally, the 
evolution and collapse of a quasi-one-dimensional vortex Airy beam in a Kerr medium are analyzed as a limiting case. 
Our study shows the possibility of controlling and manipulating the collapse, especially the position of collapse in the 
vortex Airy beam, by choosing the initial powers, vortex orders or modulation parameters. Analogous to optics, the 
results can be extended to manage the collapse of other nonlinear waves such as fluid and matter waves. 

The moments approach analysis. The propagation of a light beam in a Kerr medium is described in the paraxial 
approximation by the following nonlinear Schrodinger (NLS) equation: 

ViE-2ik^ + ^-^\EfE = 0, (1) 
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where k is the linear wave number, x and y are the transverse coordinates, z is the longitudinal coordinate, no is 
the linear refraction index of the medium, n2 is the third order nonlinear coefficient. Due to the complexity of the 
evolution of a vortex Airy beams in a Kerr medium, we apply the approach of moments [21] to study the nonlinear 



dynamics by analyzing the evolution of several integral quantities derived from the NLS. These quantities are defined 
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These quantities are associated with the beam power /i, beam width /2, momentum /s, and Hamiltonian I4; and 
satisfy a closed set of coupled ordinary differential equations [21]: dli{z)/dz = 0, dl2{z)/dz = hiz), dls{z)/dz = 
4/4(2:), dl4{z)/dz = 0, and the important invariant under evolution, Q = 2/4/2 — -^1/4, Therefore, the following 
Ermakov-Pinney equation describing the dynamics of the scaled beam width can be obtained: 
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For a vortex Airy beam, an initial field distribution can be described by [19]: 

E[x, y;z = 0) = AoAi{x/xo) exp{axx/xo)Ai{y/xo) exp{ayy/xo){x + iy)"^, 



(4) 



where Aq is the amplitude of the complex amplitude E{x,y,z = 0), is an arbitrary transverse scale. Exponential 
factors ax and ay are positive in order to ensure containment of the infinite Airy tail in the —x and —y directions, 
respectively. The azimuthal index m represents the order of vortex. The general solution to Eq. (3) with the vortex 
Airy beam as an initial field distribution can be given as 



I^{z) = hiz = 0) + 
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Eq. (5) describes the variation of the scaled beam width of the vortex Airy beam in a Kerr medium. When Q = 0, the 
rms beam width remains constant as recognized from Eq. (5) and the critical value, /f^, is derived from 2/4/2-/3/4 = 
and Eq. (2). The corresponding power can be found through Per = [21]. The expression for Per is not 

presented here due to its lengthy expression. 

Figure 1 shows the denary logarithmic vertical scale of the ratio of the critical power of a vortex Airy beam (with 
m = 1, 2) to the critical power P^ = ttcsotiq/ (712 A:^) of a Gaussian beam , where c is the speed of light in vacuum and 
^0 is the permittivity of free space. In Fig. 1 one sees that the value of the critical power of the vortex Airy beam 
depends mainly on the beam profile of the transverse distribution (besides the nonlinear parameters of the medium) . 




FIG. 1: (Color online) The critical powers of vortex Airy beams for different ax and ay. (a) m = 1; (b) m = 2. 

The critical power of the vortex Airy beam increases as the the order of vortex modes increases, but decreases as the 
modulation parameters ax and ay increase. If the initial power exceeds the critical power Per, the rms beam width 
goes to zero at a finite propagation distance, as predicted by the approach of moments. Obviously, the critical power 
is the upper bound for the collapse power [22, 23]. 

Evolution and collapse of vortex Airy beams in Kerr media. Let us examine the evolution of a vortex Airy beam 
in a Kerr medium. In the numerical calculation, we take A = 0.53/im, xq = 100/im and zq = kx'^/2 = 6cm. The 
intensity distributions of the 1st and 2nd vortex Airy beams with Pin = bP^ and Pin = 15 P^ in the focusing 
nonlinear medium at different propagation distances are shown in Fig. 2, where hereafter the intensity distribution 
is normalized with respect to its initial peak intensity. One sees that the beams with different initial powers in the 
Kerr medium propagate along the same accelerating curve trajectory as an Airy beam in free space, as expected. The 
collapse will not occur during the propagation for these cases but the evolution of intensity distributions of the beam, 
however, depends sensitively on the order of vortex and the initial power. 

As the initial power increases, more energy accumulates at various locations of the beam, leading to partial collapse 
of the beam. If the initial power is Pin = ISP^,^, the 1st order vortex Airy beam collapses at the outermost main 
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FIG. 2: (Color online) (Color online) The intensity distribution of vortex Airy beam (ax = ay = 0.1) at different propagation 
distances with initial powers in focusing Kerr medium. Upper row:the m — 1, Pin = 5Pc^ = 0.0023Pcr; and lower: m = 2, 
P^^ = 15P^^ = 0.0024Pcr. 




FIG. 3: (Color online) The intensity distribution of the vortex Airy beam {ax = ay = 0.1) in focusing medium, (a) m = l,Pin = 
15P^^ 0.0069Pcr, z 5zo; (b)m 1, Pin ISP^? 0.0069Pcr, z 6zo; (c) m 2,Pin SOP^? 0.0079Pcr, ^ 7zo; (d) 
m 2,Pin 50Pc? 0.0079Pcr, z 7.3zo; 

lobe of the beam, as illustrated in Fig. 3 (upper). No collapse of vortex Airy beam with the order of m = 2 is 
expected, if Pin = 15P^ (see Fig. 2 (lower)). The 2nd order vortex Airy beam, however, is expected to collapse at 
Pin = 50P^^, as shown in Fig. 3(c) and 3(d). In this case, the partial collapse will occur separately in two lobes in 
the beam (but, not in the outermost lobes). Numerical simulations indicate that the actual collapse power for the 1st 
order vortex Airy beam with a^^ = = 0.1 is Pin = 0.0063Pcr = ^^P^ and that for the 2nd order vortex Airy beam 
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FIG. 4: (Color online) The intensity distribution of the vortex Airy beam {ax = ay = 0.1) in focusing medium, (a) m = l,Pin = 
50P^r = 0.023Pcr, z = 0.7zo; (b) m = l,Pin = 50Pc^ = 0.023Pcr,^ = O.Szo- (c) m = l,P^n = 2500Pc? = l.lSPcr, ^ = 0.02^o; (d) 
m 1,P^^ 2500Pc^ l.lSPcr,^ 0.03^0 ; 



■ free space 
Pin=5P: 




(a) 




(b) 



FIG. 5: (Color online) The peak intensity of vortex Airy beam = ciy = 0.1 as a function of the propagation distance with 
different initial powers (a) m = 1; (b) m — 2. 
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with ax = ay = 0.1 is Pin = 0.0058Pcr = 37P^. It is interesting to see that the 1st order vortex Airy beam collapses 
at the outermost lobes, if Pin = 50P^ (see Fig. 4(a) and (b)). A stronger initial power increases the number of 
lobes that collapse at the outermost region, and decreases the propagation distances of beam (e.g. Fig. 4(c) and 4(d) 
with Pin = l.lbPcr = 2500P^). In Fig. 5, we show the normalized peak intensities as a function of the propagation 
distance with different initial powers. Although the approach of moments predicts that the rms beam width will be 
broadened when Pin < Pcr^ the results from the variation of peak intensities suggest a deformation and redistribution 
of the vortex Airy beams during the propagation. The intensity at some lobes will dominate and eventually collapse 
while the rms beam width increases or remains constant, such as for the cases Pin = ISP^? and SOP^^ for the 1st order 
vortex Airy beam and Pin = SOP^^ for the 2nd order vortex Airy beam. The results indicate that a higher initial 
power leads a shorter propagation distance before collapse. When the initial power begin to exceed the collapse power, 
the flow of transverse energy of Airy beam and vortex leads to a collapse occurring at the outermost main lobe of the 
beam. It is interesting to see that the energy accumulates to the position where the field is almost zero and finally 
a collapse occurs. On the other hand, no collapse occurs at a position where initial field mainly exists. In all the 
studies reported in the literature (for either vortex or non- vortex beams), the field never collapses at a location where 
the initial field is almost zero [3,22,23]. In the present study, however, a unusual collapse is observed for the present 
vortex Airy beam. The vortex Airy beam may collapse at a location where the initial field is nearly zero, e.g., the 
collapse point in Fig. 3(b) (cf. the zero initial field at this point, marked as point A on the first subfigure of Fig. 2). 
During the propagation, the energy accumulates at this collapse position. The higher the initial power value, the 
closer the collapse position to the beam's center and the shorter propagation distance before collapse. These results, 
therefore, provide useful information on how to spatially manipulate a collapse in experiment. If the initial power 
increases further to a certain threshold value, the beam collapses at many side lobes simultaneously, but not in the 
outermost lobes, see e.g. the collapse position in Fig. 3(d) (cf. the initial field at these points around point A in Fig. 
2 with 2: = 0). In this case, we see that the position of collapse, the number of collapse lobes, and the propagation 
distance (before collapse) also depend on the initial power. As the initial power reaches a considerably high value, 
the collapse will appear at many outermost lobes with a very short propagation distance with a higher initial power 
and a larger number of collapsed lobes. Even if the approach of moments predicts that the rms beam width will go 
to be zero at some propagation distance when the initial power exceeds the critical power Per-, the vortex Airy beam 
collapses separately at the outermost lobes as shown in Fig. 4(c) and 4(d). 

Quasi- one- dimensional vortex Airy beam. We now investigate the effect of modulation parameters on the nonlinear 



8 




FIG. 6: (Color online) The intensity distribution of a vortex Airy beam {ax = 1.5, ay = 0.05) at different propagation distances 
with initial powers (a)Pin = 1.5Pc? = 0.00043Pcr; (b) Pin = 5.5Pc? = 0.00157Pcr. The thin red line in the first plot indicates 
the position of the longitudinal cross-section. 

evolution of the vortex Airy beam. When ax = 1.5 and ay = 0.05, the beam becomes a quasi-one-dimensional vortex 
Airy beam, as shown in Fig. 6. Under a low initial power of Pin = l.bP^, the beam propagates along the accelerating 
trajectory with a feature similar to the propagation behavior of an Airy beam in free space, the field distribution is 
evidently different with that of the beam in free space, and the beam does not collapse as shown in Fig. 6(a). By 
increasing the initial power from Pin = 1.5P^ to 5.5P^^ the off-center part of the beam partially collapses at certain 
propagation distance although the propagation trajectory is the same as that of Pin = l.bP^, (compared with the 
zero initial field of this collapse point to the point A on first subfigure in Fig. 6(b)). On the other hand, no collapse 
occurs at the position where the field originally exists as shown in Fig. 6(b). The underlying physics are similar to 
the collapse behavior of a vortex Airy beam. 

In conclusion, we have systematically studied the collapse dynamics of vortex Airy optical fields in a Kerr medium. 
The collapse of a vortex Airy beam requires a higher energy than that of a non- vortex beam. Although vortex Airy 
beams with different initial powers propagate along a similar acceleration curve trajectory in the Kerr medium as 
a non- vortex Airy beam propagates in free space, the evolution and appearance of the collapse in a vortex Airy 
beam exhibits unusual features. The collapse can occur in the position where the initial field is almost zero. Our 
study shows that the collapse can occur at different locations other than the central point, and the locations can be 
controlled through choosing appropriate initial powers, the order of vortex, and the modulation parameters. This 
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